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Abstract 

We investigate autonomous perturbations on the orbits of LISA, namely the effects produced 

by gravitational fields that can be expressed only in terms of the position, but not of time in 

i^ the Hill frame. This first step in the study of the LISA orbits has been the subject of recent 

papers which implement analytical techniques based on a "post-epicyclic" approximation in 

br^ the Hill frame to find optimal unperturbed orbits. The natural step forward is to analyze the 

r/^ perturbations to purely Keplerian orbits. In the present work a particular emphasis is put on 

the tidal field of the Earth assumed to be stationary in the Hill frame. Other relevant classes of 

autonomous perturbations are those given by the corrections to the Solar field responsible for a 

slow precession and a global stationary field, associated to sources like the interplanetary dust 

I""! or a local dark matter component. The inclusion of simple linear contributions in the expansion 

^^ of these fields produces secular solutions that can be compared with the measurements and 

I possibly used to evaluate some morphological property of the perturbing components. 

1 Introduction 

> 

\Q Recent papers have revived the analytical study of the orbits of the three spacecrafts of the LISA 

t*^ 'constellation' [H [21 13]. Although the topic has been explored and analyzed in some detail with 

^0 the usual numerical methods and many perturbing effects were already considered and dealt with 

in E] 5 we think it can be useful to pursue an analytical approach, that can be of help in gaining 
f^ insight on the relevance and hierarchy of the various perturbing effects. 

C^ The main results obtained so far include the possibility of realizing first-order stable solutions 

'^ of arbitrary shape in both 'standard' planes at angles near ztGC with respect to the ecliptic and the 

>> optimization of the configuration which produces a peak-to-peak variation in arm length of about 

k> 40000 km and a maximum rate of change of arm length of about 4 m/s. These figures represent the 

?H least possible amount of 'flexing' of the arms of the nominal equilateral triangular configuration and 

refer to a purely Keplerian motion in the field of the Sun only, on suitably phase-shifted Earth-like 

orbits. They can be obtained either by a trial/error procedure based on the direct computation 

of the exact Keplerian orbits [3] or by a perturbative analytical study. It can be shown [2] that, 

due to the small eccentricity and inclination of the orbits, a series expansion up to second order 

in a suitable perturbation parameter differs from the exact Keplerian solution by less than 0.03% 
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making the method of analytical series expansion a viable route to study more complex aspects of 
the motion of LISA. 

There are in fact several sources of perturbations affecting the motion of the three spacecrafts on 
the unperturbed elliptic orbits. It is easy to convince ourselves that all of them are weak enough to 
produce very small effects over the time of the mission. However, in view of the need of an accurate 
knowledge of distances and Doppler shifts in the implementation of the Time Delay Interferometry, 
knowledge, analysis and, possibly, mitigation of all perturbative effects are a necessary task. A 
preliminary study can be profitably done with those analytical techniques mentioned above, in 
parallel with more accurate but painstaking numerical simulations [H [5] . 

The most useful setting in which to perform the analytical study is that based on the Hill system, 
recalled in sect. 3. This is a frame rotating with a constant angular velocity set by the mean motion 
along a reference orbit. The dynamics in this frame are determined by a tidal approximation of 
the gravitational potential of the Sun and other sources. The far tide approximation [8] consists 
in taking only the first order term in the field expansion around the origin of the frame. In order 
to be consistent with the above mentioned series expansion of the Keplerian orbits, it is sufficient 
to truncate the expansion to the second ('octupolar') order of the field of the Sun. The same 
procedure can then be applied to the Earth, planets, etc. 

These perturbations can be classified as autonomous and non- autonomous. In the first class fall 
the fields that can be expressed only in terms of the position (and, possibly, of velocity), but not 
of time in the Hill frame. In this paper we shall focus on the autonomous perturbation effects on 
LISA orbits. The most relevant example in this class is the field of the Earth, that in the Hill frame 
is practically at rest. The equations of motion can be cast in the form of a system of oscillators 
subject to a linear vector potential (due to the Coriolis acceleration) and coupled through higher- 
order terms. The linear system can be solved exactly, the nonlinear couplings are then included 
perturbatively. In sect. 4 we investigate the configuration, i.e. the phase of each LISA spacecraft 
orbit, at the closest approach to the Earth, in order to see the effects on all the flexing and Doppler 
shift indicators, in particular when the time base-line of the mission is longer than usually assumed. 

The effects of global stationary fields can also be studied by their series expansion in the Hill 
frame. The linear part already provides useful information. A first interesting example, that we 
address in sect. 5 is given by the corrections to the gravitational potential of the Sun, treated as a 
classical spherical mass. The general relativistic correction in the weak field limit, responsible for 
the classical relativistic perihelion advance, can be obtained by adding a term scaling as the inverse 
cube of the distance and parametrized by the Schwarzschild radius. Another interesting correction 
term is provided by the intrinsic quadrupole moment of the Sun. With the current values of the 
parameters, the effects produced by both these terms are probably within the detection capabilities 
of LISA. In sect. 6 we discuss the tiny effects of the interplanetary dust or a local dark matter 
component that can in principle be traced in the spectrum of the relative motion between the 
spacecrafts and can be used to evaluate some morphological property of the perturbing component. 

In general, the perturbative effects due to all other bodies in the Solar System have to be 
expressed by fields explicitly dependent on time and therefore fall in the non-autonomous class. In 
sect. 7 we conclude with some perspective on the study of this second class of perturbations. 



2 Unperturbed Keplerian orbits 

We briefly recall the main results, obtained by other authors, that are useful for our work. We 
begin by writing the unperturbed motion of the three spacecrafts SCk, k = 1,2,3, on Keplerian 
orbits in the field of the Sun alone. In solar barycentric coordinates, let us call Xk,yk,Zk the 
positions of the SCk along the three orbits. They are given by [1]: 

2tt 2tt 

Xk = Xi(V'fc)cosy(A:-l)-yi(V^fc)sin— (fc-1), 

27r 27r 

n = Xi(V^,,)sin— (fe-l)+yi(V^fc)cos— (fe-l), (1) 

Zk = Zi(^pk), 

where t/^k is the eccentric anomaly of SCk (see below, Eqj6|. We then introduce the natural series 
expansion parameter 

i I 

a = — ~ — , (2) 

2R 60' ^ ^ 

where ^ = 5 x 10^ km is the nominal length of an arm (the side of the LISA triangle) and R = 
1.5 X 10*^ km = lUA is the semi-major axis of the Earth orbit. The authors of ref [Ij choose to fix 
the eccentricity e (the same for all SC's orbits) 



/ 4 4 

e = i 1 + -a"^ -\ — ^a cos — 1 (3) 

by the choice of a and 6, the inclination of the plane of triangle with respect to the ecliptic, 

e = ±'^ + 6 = ±'^+a6i, (4) 

where 6 (alternatively 5i) denotes the additional tilt angle with respect to the canonical value ib|. 
Finally, the inclination i of the orbital planes with respect to the ecliptic is given by 

2 a sin 9 

tarn = —pz K . 5 

V3 1 + ^acos6l 

We can now write the orbit of one 'reference' spacecraft, e.g. SCi, in its parametric form: 

Xi = R{costl'i + e) cosi, 
Yi = i?v 1 — e^ sin^i, 
Zi = R{cosipi + e)smi, 

The orbits of SC2 ed SC3, are obtained starting from that of SCi and performing two rotations 
around the Z axis, of an angle 27r/3 and 47r/3 respectively, leading to eqs.Q. The eccentric anomaly 
as a function of time is given by the solution of the Kepler equation. Actually, we now have three 
anomalies, given as solutions of the three equations 

2tt 
ipk + esini/jk = ^t- —{k -I) = cl)k, k = 1,2,3 (6) 



where the mean motion is ri = 2tt/T = \/ —jw^- If we measure time in years, then T = 1 and 
SI = 27r. It is also useful to introduce the phase shifts 

27r 
ak = —{k-l), (7) 

so that (f>k = Vtt — <Jk. The distances among spacecrafts are clearly given by the three lengths 



Ljk = ^{X, - Xkf + {Y.J - n)2 + {Zj - Zfc)2, j, k = l,2, 3. (8) 

These distances are not constant in time, causing what is usually referred to as 'flexing', even for 
purely Keplerian orbits. A simple way to see this, is to consider that the instantaneous distance 
from the Sun are different from one spacecraft to the other. The flexing is quite sensitive to the 
inclination of the plane of the triangle. The value 5i = 5/8, giving the angle ~ 60° 35' 45" is 
proved to be optimal by Rajesh Nayak et al. j2]. The maximum flexing of the optimal configuration 
is about 40000 km and the maximum rate of change of arm lengths is ±4 m/s. The reduction 
of the relative velocity and related Doppler phase shift with respect to the case considered in [2] 
((5i = 0) amounts to a factor 5.5. However, it should be noted [6] that flexing figures comparable 
to the optimal ones can be obtained simply with slightly different choices of the Keplerian orbital 
parameters [H [H] . The large amount of flexing considered in ^ with (5i = is determined by the 
choice implicit in definitions ([3]-[5]). 

A nice way to appreciate the improvement provided by the optimal choice of the tilt angle is 
that of comparing the spectral amplitudes of the differential relative motion. Following ref.[TT], we 
can compute the spectra of the difference 

ALfcj = Lij - Lik, i,j, k = 1,2, 3. (9) 

In FigjT] the two curves represent the modulus of the Fourier transform of AL32 tapered with a 
Blackman window. The spectrum is computed on a time span of 30 years chosen just to get a 
good resolution in frequency. The optimal curve has the second harmonic (frequency 2yr~^) seven 
orders of magnitudes below that of the standard solution with 5i = 0. The higher harmonics are 
much smaller or even unnoticeable. In both solutions the harmonics 3, 6, 9, ... are absent in view of 
the triangular symmetry of the configuration. 

3 The Hill system 

A perturbative analytical approach is convenient when we have to analyze in the simplest way the 
geometry of orbits and the most relevant perturbations. The Hill reference system (see, e.g. [3 [9], 
also called as Clohessy- Wiltshire system [14] ) allows us to reduce the equations of motion to those 
of a set of coupled oscillators. The Hill equations properly said consist of the linear part of this 
system which is readily solved. Non-linear and non-autonomous terms can be added and solved 
perturbatively. In the specific instance, the system is at rest with the barycenter of the constellation 
and so an observer sees the Earth at rest in it. The new reference is introduced according to the 
following steps: 

1. Start by considering a circular orbit of radius R around the Sun and pass from the inertial 
barycentric system to a new system with its origin rotating with the reference orbit and with axes 
X = xi,y = X2, z = X3. xi is directed radially opposed to the center, X2 is in the direction tangent 
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Figure 1: Logarithm of the modulus of the Fourier transform of the difference L12 — L13: standard 
tilt {9 = tt/3, di = 0) (continuous line); 5i = 5/8 (dashed line). 

to the motion and X3 is perpendicular to their plane. The system is non-inertial and the equations 
of motion are 



X = -V$ -2Q,Ax-nA{nA{x + R)). 

$ is the total gravitational potential and fl = (0,0,^1.), R= (R, 0, 0). 
2. The potential is approximated according to the expansion 



(10) 
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x^xx + 0{\xf), 



(11) 



a;=0 



x=0 ^ ^X^UX\ 

where the sum is implicit over the repeated indexes A, /i = 1, 2, 3. The "far tide" (or Hill) approxi- 
mation [8] consists of retaining only the first order term into the expansion of the gradient of the 
potential. Taking 
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the Solar far tide approximation gives 



GMp 







V(x + RY + y'^ + z'^ 



5$ 
dx 



x=0 



GMq 2 9$ 



x=0 



dz 



(12) 



(13) 



x=0 



and 



dx'^ 



2GMp 







x=0 



R3 



-2Q^ 



92$ 

(9y2 

5 



x=0 



92$ 

aT 



GMp 







x=0 



i?3 



n' 



(14) 



all others zero. The next order of approximation {octupole) is: 



^3$ 






qgmq m^ 



dxdy'^ 



x=0 



x=0 






dxdz 



all others zero. By using these in (10) we get 
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(17) 



These are the equations for a second order approximation of the purely Keplerian motion. Any other 
perturbing force can be derived by an analogous series expansion of the corresponding potential and 
added to this system. Each component, coherently with a perturbative approach, will be a power 
series in the coordinates with suitably 'small' coefficients s = {ex,£y,£z)- Recalling the definition 
of the parameter a in ([2| and introducing explicitly only additional autonomous perturbations 
represented by the perturbation vector f{x,e), the system (10) at second order can be written in 
the form 



X - 2^y-3n^x+^^{2x^-y^-z'^)-Ux,e) = 0, 

y + 2ite —xy-fy{x,e) = 0, 

2 Q^^a „ , X „ 

z + Hz — xz — jz[x,e) = 0. 



(18) 



Since a is small, the nonlinear terms in the multipole expansion of the Solar field can be included 
in the perturbations. This is an important remark, since in general the presence of quadratic or 



higher-order terms in the system ( 18 ) makes it non-integrable. The most effective approach is that 



of finding the exact solution of the linear part and using it recursively to find approximate solution 
for the full non-linear system. 

The perturbation vector f{x, s) could even start simply with linear terms of the form ex + . . . 
giving explicit small changes to the unperturbed zero-order solution. However, it is more natural 
to consider them of magnitude comparable to terms of order a^ or higher. These and possibly 
higher-order perturbing terms propagate in the full solutions through the non-linear coupling with 
the other components of the field expansion. The solution of the system can therefore be written 
in the form of a series expansion 



.(0) 



.(1) 



(2) 



+ ax) + ExXl + 



Xkit) 



(19) 



where the value of the upper index denotes the order and with the index k = 1,2,3 we denote as 
usual the three spacecrafts. 



The 'zero' order solution is obtained by putting a = e = in (18). It amounts to be 



xi{t) = 2(^2xi(0) + ^yi(0)j +^xi(0)sinm- (^3xi(0) + ^yi(0)jcosfit, 

yi(t) = yi(0)-^ii(0) (20) 

-3 (2Qxi(0) + yi(0)) t + -xi(0) cos Vtt + 2 ( 3xi(0) + -yi(0) J sin VLt, 

Z]^{t) = zi(0)cosf]t + -^sinOt. 

The solutions for the other two spacecrafts are obtained by substituting Vtt with 0^ as in definition 
(|6]). However, the solution actually employed in the perturbation expansion (the epicydic motion) 
is obtained by choosing initial conditions in order to remove constant and secular terms [1]: 

(0) ^ A (0) ^ ■ A (0) ^ A ^o1^ 

4 =^cos(/>fc, yl' = --^sm(t)k, z^. ' = -coscpk- (21) 

It corresponds to harmonic oscillations around the 'reference' x = circular orbit. 

The order 'one' (in a) solution with the integration constants chosen in order to comply with 
the LISA constraints is [2] 



„(i) 



5^ ^ /^l . 1 



,(1) ^ / 1 



C = -^^-Sijsm(Pk + - sin 2^k, (22) 






^k^ = ^-r^(l-'5i)cos</.fc--^cos2</.fc. 



The solution 



Xk{t) = xf^ + ax^^^ (23) 

with these explicit forms is a very good approximation of the exact Keplerian orbits, with an error 
less than 0.03% on the prediction of the flexing of the arms [21 13]. In other terms, the plots of 
Fig jl] can be reproduced using solution (23) with a remarkable accuracy. This result proves that 



this solution is a reliable starting point for the study of additional sources of gravitational distur- 



bances. From the analytical point of view it is worth mentioning that also system (18) with just 



f{x, e) = (of which (23 ) is a first order approximate solution) is not integrable unlike the original 
Kepler system which is clearly integrable. Is a common lesson of modern analytical mechanics 
that truncated solutions of non-integrable approximations can be more useful than formally exact 
solutions of their integrable counterparts. 



4 The Earth effect 

As mentioned above the Earth is practicahy at rest in the Hill system. The Earth coordinates are 

xe = -R{1 - cos 20°) ~ -9.0 x 10*^ km, y© = i? sin 20° ~ 5.13 x lO'^ km, z© = 0. 

We can therefore consider its effect as a constant + linear perturbation to add to the equations 
of motion. If d© = \/x'^ + y^ ~ 5.2 x lO'^ km, is the distance from the origin, we can add the 
perturbative parameter 



where we used the ratio 
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(24) 



The equations of motion can now be written in the form (18) with the perturbation vector f{x, e) 
given by 



Proceeding as assumed in eqs.(19), we may look for a solution that can be linearly superposed to 



(23) so that we have 

Xk{t) = xf + a4^^ + e^x^^K (26) 

The explicit solution is [7]: 



„(2) 



Ah 5i 

2— +xq + 2y©(/)fc + Bk cos (l)k +Cksin(t>k - —i=(l)ksm(j)k, 



yf = -(3^ + 2x©)0fc-^y©<Ai+ (27) 

5£ \/3£ 

hk cos (/)fc + -—- sin 0^ + 2(Cfe cos (/)fc - Bk sin 0^) + Dk, 



2x/3'^ 2 

(2) ^ 

4 = Ek cos (f)k + Fk sin (f)k - j4>k sin (f)k, 

where the Ak,..-,Fk are constants determining the initial conditions. Recalling that the mean 
anomalies <J3k are just phase shifted rescaling of time, in these solutions we see the presence of 
additional secular terms. In particular, the y coordinate has a term quadratic in time that cannot 
be eliminated by any choice of the initial conditions. 

Dhurandhar et al.[3 propose to inject LISA on orbits with initial conditions devised to reduce 
to a minimum the Earth effect. In view of the secular terms, it seems a good strategy that of 
minimizing the disturbance at 'mid mission': they assume a mission duration of three years, so 
that this 'optimal time' is at 1.5 years after starting the operation activity. The authors of rcf.[7j 
introduce in the mean anomalies (pk a phase shift to accommodate for different starting times and 




Figure 2: Distances Lij vs time including the Earth influence for a period of 6 years: L12, continuous 
line; L23, long dashes; L13 short dashes. 
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Figure 3: Rates of change Vij = Lij vs time including the Earth influence for a period of 6 years: 
1112 , continuous line; t;23, long dashes; V13 short dashes. 



compute the resulting effect on the chosen period. However, the solution is invariant under time 



translation in view of the autonomous nature of the system, therefore, expressions (27) can be 
profitably used at arbitrary initial time. If this, for simplicity, is put equal to zero, the values of 



the constants appearing in (27) are 
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Ck 



D, 



Eh 






13^ 

8\/3 



, cos ak - 27/0 sin ak 



V3i 



, sin ak - 2y© cos ak + £ 



3 sin 2c7fc 



- COs2c7fc, 
lOcTfc 



8^/3 



(28) 



A 



2x©(Tfc 



v3£(7fc cos ak H — ^ sm ak , 






SVC? ak, 



{2ak -sm2ak) 



where the ak are the phase shifts introduced in (JTl). 

Experience gathered on many space experiments of the past allows some optimism on expecting 
the LISA mission to last longer than the nominal duration of three years. A longer lifetime also 
allows us to better investigate the secular (non periodic) behavior of the orbits. In order to clarify 
the effect of the Earth tidal field, we have therefore chosen to consider an "optimistic duration" 
and to double the expected mission time: in Figj2] the three arm lengths time are plotted over a 
period of 6 years with the solutions (26) and the initial conditions above. In Figjslthey are plotted, 
on the same six years period, the relative velocities Vjk = Ljk- Now we can see that, if we limit 
the observation to the three years centered around zero, both the flexing and the relative velocity 
have an almost constant envelope. However, in the long run the envelope tends to linearly grow for 
both quantities and is therefore appropriate a comparison with the 'optimal' solution illustrated in 
section 2. In a three- year (±1.5) mission we have a flexing of 60000 km which increases to 80000 km 
in six (±3) years: on this longer time span we arrive at a factor 2 of increase of the flexing with 
respect to the purely Keplerian solution. The maximum rate of change of arm lengths is ±5.5 m/s 
in the ±1.5 years time span but increases to ±10 m/s in the ±3 years time span, a factor 2.5 of 
increase with respect to the Keplerian solution. 

This behavior is due to the global dynamics in the coupled fields of Sun and Earth as described 



in the rotating frame. As can be seen from the solution (27) with the constants determined by (28) 



the secular terms always produce a drift away from Earth: this is linear in x and, as remarked 
above, quadratic in y. The secular drift in x of the three SCk amounts to ±70000 km on the ±1.5 
years time span, whereas the motion in y is of initial approach, a state of rest and a subsequent 
departure: on the same time span, the drift in y is of 5 x 10^ km. So, the best that we can do is a 
choice of initial conditions such that LISA is a little further out at start, approaches the Earth, is 
at rest ad mid mission and departs after that. The physical reason of this behavior is evident when 
we consider the nature of the effective potential. LISA is well outside the Roche lobe associated 
to the Earth; therefore it cannot 'fall' (as a stone or a projectile [7]) towards it. Rather, it moves 
on a slow merry-go-round and feels a small Coriolis apparent force. The combined parabolic-drift 
motions have therefore quite a relevant role in increasing the flexing on a certain time base. 

However, the possibility exists of an appreciable reduction of the flexing due to the Earth 
tidal field over a limited time span by a suitable shift in the relative phase between the epicyclic 
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Figure 4: Distances Lij time including the Earth influence for a period of 6 years with an extra 
7 = — 7r/2 phase: L12, continuous line; L23, long dashes; L13 short dashes. 



motion and the linear perturbation. The scenario depicted above can in fact be changed if the 
constellation gets to the point of minimum distance from the Earth with a different orientation. 
This situation can be described either by modifying the phase shift in the initial conditions of the 
linear perturbation or by putting this extra shift directly in the definition ([7|. Let us modify this 



m 



2tt 



<7fc 



ik-l)-jk, k = 1,2,3 



(29) 



using the angles 'jk to change the orientation of the triangle at a given time. We limit the analysis 
to identical angles for all the SC^, however a non-rigid rotation could also be studied. The standard 
zero order solution used above, with 7 = 0, of ( 21 ) corresponds to a configuration with the spacecraft 
1 at the 'highest' point {zi = i/2) over the ecliptic at time zero (let us say at mid mission). The 
spacecraft 1 being at zi = with either yi = +£/\/3 or yi = —£/y/3 correspond to the two extremes 
in which it is closest or farthest from the Earth at time zero. It is interesting to examine these 
two situations compared with that above: they are given by putting 7 = ±7r/2 in (29). In FigJ4] 
the first case (7 = — vr/2) is illustrated by plotting the three arm lengths time over the period of 6 
years: in this case in which SCi is the one closest to the Earth at time zero, the two arms 12 and 
13 exhibit a variation not larger than the optimal one (40000 km) for the first half-mission period 
of 3 years. It can be shown that in the second case {SCi farthest from the Earth at time zero) the 
plot is the time reversal of that in Figj4j so that the variation of L12, -L13 is kept to a minimum in 
the second 3 years. The price for this reduction is that the flexing of the third (23) arm is higher, 
but only slightly with respect to the level found with 7 = (80000 km) on the same time span. 

Another way to examine the effect of the Earth is by comparing the differential relative motion 
computing the difference AL32 = L12 — L13 and the corresponding spectral amplitudes. In FigJ5] 
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Figure 5: Spectrum of the difference AL32 = L12 — -L13 for the case 7 = 0. The dashed hne is the 
Keplerian configuration and the continuous Une represents the addition of the Earth perturbation. 

the modulus of the Fourier transform of AL32 is shown together with the optimal Keplerian case 
of Figll] (no Earth contribution). The continuous line includes the contribution of the Earth in 
the case of 7 = 0. As before, to get a good resolution in frequency, the spectrum is computed 
on a longer time span (30 years). The effect of the Earth is evident from the appearance of the 
harmonics 3, 6, 9, ... due to the fact that the triangular symmetry of the configuration is now broken. 
The configurations with the phases 7 = ib7r/2 display the same spectral lines with slightly lower 
values at higher frequencies. 

The analysis of these and other possible configurations could be of a certain practical relevance 
should it be necessary to privilege the proper working of only a pair of arms [16] . 

5 The relativistic precession and the effect of the Solar quadrupole 

In the list of small perturbative effects given by fixed fields naturally fall those connected with 
the not exactly Keplerian interaction with the Sun. Two important examples are the relativistic 
correction of the classical field and the perturbation due to a possible quadrupole moment of the 
Solar gravitational field. The first produces the slow precession responsible for the advance of the 
perihelion; the second produces a more complex variation of orbital elements but of (supposedly) 
much smaller amount. 

It is tempting to see if LISA is capable to detect small effects due to the violation of the 
spherical symmetry of the Solar field, such as we have assumed so far. The small eccentricity and 
the distance from the Sun of Earth-like orbits makes the amount of the effects quite small, but 
worth investigating in view of the high sensitivity of LISA to small differential variations of the 
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arm lengths. On this ground there is the recent proposal by Polnarev et al.[15] of detecting with 
LISA the low frequency g-mode Solar oscillations. 

The relativistic correction can be mimicked by the presence of a small additional term in the 
potential of the form 



where 

GMq 

is the Schwarzschild radius of the Sun and L is the angular momentum of the orbit related to the 
other parameters by L^ = GMqR{1 — e^). The Solar quadrupole is associated with a small term 
of the form 

^Q = -3 2r^2 ' (30) 

where J2 = Q x 10"''^ characterizes the mass quadrupole of the Sun, Q being a parameter of order 
1 according to current estimates [18]. Rq ~ 7 x 10^ km is the radius of the Sun and the radial 
distance is expressed in terms of the coordinates as 

r = \JX'^ + Y'^ + Z'^ = ^J{R + xY + y'^ + z'^. 

An approximate evaluation of the two effects can be attempted in the Hill frame with the usual 
hypothesis of small perturbations of almost circular orbits. The potentials with the perturbing 
terms introduced above are expanded in series about the origin of the Hill frame and the first 
and second order terms are retained and added to the Keplerian terms. We remark that a small 
constant term appears in the x component because the centrifugal force in the rotating frame 
exactly balances only the Keplerian monopole ak, r = R: the extra term is responsible for the 
classical precession effect common to all spacecrafts. The linear terms are given by the different 
quadrupole terms arising in the expansions: these give significant contributions in the form of 
differential effects. The equations of motion are of the same form for both effects and can be 



written in the form (18) with the perturbation vector f{x,e) given by 

f, = -e,n''iR-4x), fy = -eyn^y, f, = -e,n''z. (31) 

The Schwarzschild correction is obtained if the coupling constants are given by 

ex = ey = e^ = £s = 3^ ~ 3 X 10-s (32) 

whereas the Solar quadrupole is included when the coupling constants are given by 

1 3 ^ Z' Rg\ ^.00^.. in-12 



ey=3^. =^Q=2"^Mlf ) -3.3QxlO-^^ (33) 



The structure of the system is the same as in the case of the Earth (cfr. Eqs.(18) with the 



perturbation vector given by (25)) and the strength of the perturbing terms can be evaluated 
by comparing the above constants to £9 ~ 7.3 x 10"^. For the consistency of the approach we 
preliminarily check the prediction concerning the overall precession. The constant force appearing 
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J, 



5^ 




/(years) 



Figure 6: Difference AL^^*^ — AL •^'^ time for a period of 5 years as a consequence of the differential 
effect due to the relativistic precession: the continuous hne is for jk = 23, the long dashed line for 
jk = 13 and the short dashed line for jk = 12. 




?(years) 

Figure 7: Difference ALgg*^ — AL32^ in units of Q time for a period of 5 years as a consequence 
of the differential effect due to the Solar quadrupole: the continuous line is for jk = 23, the long 
dashed line for jk = 13 and the short dashed line for jk = 12. 
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in the x component is a factor R/i = (2a) ^ = 30 stronger than that hnearly depending on the 
coordinates. Therefore we can assign to the solution the following structure 



Xk{t) = x'^^ + ax'^^ + ^a;(2) + x^^\s). (34) 



In the second-order term there is no index k because the solution satisfying proper initial conditions 
on the epicycle is the same for each spacecraft. The structure of this term has to comply with the 
natural condition that the radial oscillation be symmetric with respect to the circular reference 
orbit whose radius is slightly reduced by the constant component of the force. The choice of initial 
conditions 

I £ £ 

Xk{0) = —^coscTfc, yk{0) = y=smak, 2^(0) = - cos 0"^; 

n£ , , m f e^ \ , , ^£ 

Xk{0) = — ^smcTfc, yk{0) = — y=lcosak-—^j, za:(0) = — smcjfc, 

gives the correct solution 

x(2)(t) = _-(i_cosm), y^^\t) = £ (nt - - sinnt) , z'-^^t) = 0, (35) 

since the constant part in x{t), —ex£{Qa)~^ = —SxR/S, accounts for the change of radius of the 
reference orbit. An estimate of the advance of the perihelion after a revolution is given by 



2a R 



A^^^'-^Jl^=27rey. 



Using Eqs.(32) and (33) we get respectively 



^Vs — — 75~ [rad/orbit] 
R 

for the relativistic advance of the perihelion and 

2 



A(^Q ~ 37r J2 f -^ j [rad/orbit] 



for the corresponding advance due to the J2 term. These results are in agreement with standard 
findings in the limit of small eccentricities [19] making us confident about the reliability of the 
solution. The higher order terms are 
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xf ) (t) = -2xf it) = -2xf (t) = ^ (1 - cos nt) , 
yj^) (t) = -^y3eyi {nt - sin nt) , 

y?it) = syi (-^ + ^nt- COS (^nt-^^Y (36) 



y(3) (i) = Syel^+—Qt + COS ( m 



,(3),,^ _ ^ J^ ^^n,^^^Jn,_^\\ 



zf^(t) = _^e^£ ( [7tsm(f^t- ^(/c- 1) ) - ;^^'^ ^^sinJ^t 



To better show the effects produced by these perturbations we can use the quantities 

AL^f -A4f\ (37) 

where AL^j = Lij — Lik and the arm-lengths Ljk are computed using the full perturbed solution 



(34) in AL^'j^ and including only the Keplerian solutions in AL -^^ . 



In FigJ6]we show these quantities time for a period of 5 years as a consequence of the differential 
effect due to the relativistic precession: the amplitude of the arm length variations can be as large 
as 10 km. In Figj7]we show the same quantities due to the differential effect of the Solar quadrupole: 
the amplitude of the arm length variations reaches 1.5 m times Q. The modulus of the Fourier 
transform of AL^^ — AL.^ is shown in Figlsl 

In principle, a precise evaluation of all other effects (Schwarzschild precession included) would 
allow a measurement of Q (and J2) from the residuals. Indeed, Figjslshows that the relative content 
of higher harmonics with respect to the fundamental is larger for the Solar quadrupole precession 
than for the relativistic one. However, this measurement is hindered by the absolute magnitude of 
the first effect being orders of magnitude smaller than the latter. 

6 A global stationary potential 

The presence of a stationary global potential in the Solar System can be due to different kinds of 
diffuse matter distributions. Main candidates are the interplanetary dust and a local dark matter 
component |101 lllj . Our concern here is to show that this kind of gravitational source is readily 
described as an autonomous perturbation and can be analyzed with the methods illustrated above. 
We consider a very simple model in which a spheroidal distribution can be assumed with rotational 
symmetry around the (barycentric) Z-axis and mirror reflection symmetry with respect to the 
X — Y plane. Another apparently very crude assumption is that of assuming a linear gradient field 
as a second order perturbation in our system of equations of motion. This hypothesis is however 
viable since almost circular orbits are probing the extra field in a small radial range in which it is 
slowly changing and what can be of practical interest is the possibility of detecting an anisotropy in 
the field, possibly due to the ellipticity of the spheroidal distribution. The hypothesis is sustained 
by numerical computations with the interplanetary dust [TT] showing that its perturbing effect is 
weakly affected by the steepness of the radial component of the density function. 
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Figure 8: Spectrum of the difference ALgg*^ — ^L^ : the continuous hue represents the effect due 
to the relativistic precession; the dashed hne that due to the Solar quadrupole (with Q = 1). 



The system of equations that we use is of the form (18) with the perturbation vector f{x,£) 
given by 



/^. = -e-D^^x, fy = -e-D^'^y, f^ = -Pe-o^'^z. 



(35 



Here ^d is the small parameter associated to the diffuse component and we have assumed a 
spheroidal distribution parametrized by a shape parameter /3. In the very simple setting sketched 
above, we may assume 



£d 



A-kGpd 



(39) 



302 ' 

where p-Q is a suitable average over the density distribution. The solution is still of the form ( |26[ ) 
where in the second order terms we add the contribution due to the solution 



(2) 



(2) 
VkT) 



A' 5i 

2-^ + Sfc cos (pk + C'l, sin (/)fc - — /-'/'fc sin (t)k, 

A' 5i / 3 

-3-^0fc + 2(Cfc cos 0fc - B'l. sin (pk) + D'^ - — = f (f)k cos 4>k - -= sin (/)fc ) , (40) 



4d = K cos 0A: + K sin (l)k- j<t>k sin (/)fc 



and the arbitrary constants have a prime to distinguish them from those of the solution for the 
Earth. With the natural condition of putting to zero both the displacement and its time rate at 
'zero' (half-mission) time, we get 
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, 3sin2(TA; -10(Tfc , r- , 4^ ■ lA^\ 

Cfc = £ — ^= , Dfc = -\/3£(TfcCosafc + — ^smcTfc, (41) 

£'^ = --sin^o-fc, F^ = -- (2crfc -sin2(Tfc). 



The components in ( 40 ) , multiplied by ed and added to the solution ( 23 ) , provide the evolution of 
the arm- lengths under this perturbing effect. To show the small effects of the diffuse component, 
Cerdonio et al. fTT] (whose approach is purely numerical) plot the following quantities 

5L,, = ALD-4f, (42) 

where the arm-lengths LjAid are computed using the full perturbed solution and the Lj^ are those 
obtained with the optimal Keplerian solution. 

In Figj9] these relative distances are plotted in the case in which the average density is taken 
to be pd = 9.6 x 10"^'' kg jrr? as estimated by models of the meteoritic complex p!7] and a shape 
parameter /3 = \/3 (spheroidal distribution): we can compare these plots with those for the a = 
case of reference [H], corresponding to their most homogeneous distribution. Our solution is in 
good agreement with that obtained with a numerical integration of the Gauss equations [11]. In 
FigjTO] we have plotted the Fourier spectrum of two different realizations of the difference bLj,2 



of Eq.(42). In ref.[llj an analogous plot has been obtained by the discrete Fourier analysis of 
the numerical solution and we can again appreciate the good accuracy of this simple analytical 
approach. We can see how the solutions suggest the possibility of assessing the shape of the matter 
distribution through the spectral analysis of the arm lengths variation. However, we have seen in 
the previous sections that other autonomous perturbations effects leave on these observables marks 
at the same frequencies but with higher amplitude. 

7 Conclusions and future work 

We have reported on a systematic study of the autonomous perturbations of LISA orbits, namely 
those stationary with respect to the center of the LISA constellation, i.e. constant in the Hill frame. 
This is the simplest step in the study of the perturbations on the LISA orbits. 

As in other recent papers [U El [3], we implement analytical techniques based on a 'post-epicyclic' 
approximation in the Hill frame. This work is intended as a preliminary study aimed at gathering 
an overall picture and more physical insight to the issue. In fact, the analysis of changes of the 
ideal Keplerian orbits may suggest in principle possible changes in the choice of initial conditions. 
Accurate, dedicated numerical investigations [H El |9] are then required in the assessment of the 
exact amount of each effect. However, a numerical optimization depends, in the case of independent 
orbits, on eighteen parameters and is yet to be fully explored |4J. A systematic comparison of our 
analytical approach with these solutions is a further step that must be undertaken. 

A particular emphasis was put on the tidal field of the Earth, assumed to be stationary in the 
Hill frame. The solution of the equations of motion under the perturbation gives a combination 
of a reversed parabolic motion in the tangential direction and a drift in the radial direction. This 
sheds light on possible implications on the choice of initial conditions when the time base-line of the 
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Figure 9: Relative distances SLj^ time due to a spheroidal (/3 = v3) distribution of interplanetary 
dust, over a time-span of 5 years: the continuous line is for jk = 23, the long dashed line for 
jk = 13 and the short dashed line for jk = 12. 
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Figure 10: Spectrum of the relative distance 6Ljk for two values of the shape parameter f3: 13 = v3 
(continuous line) and (3 = 1 (dashed line). 
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mission is longer than that assumed in previous papers. We propose a more systematic investigation 
of the phase choice affecting the orientation of the constellation at the moment of closest approach 
to the Earth that could be of a certain practical relevance in the case it should be necessary to 
privilege the proper working of only a pair of arms. 

Another relevant class of autonomous perturbations is that given by a global stationary field, 
associated either to the classical or relativistic corrections to the spherically symmetric Solar po- 
tential or to sources like the interplanetary dust or a local dark matter component. The inclusion 
of a simple linear contribution representative of these fields produces a secular solution that can 
be used to evaluate physical parameters of the perturbing force. The relativistic precession and 
the more complex motion due to a Solar quadrupole contribution are most probably within the 
detection capabilities of LISA. More problematic, with the figures obtained, is a detection of a 
possible flatness of the perturbing spheroidal component. 

The natural step further on this research line is to attempt the analysis of non-autonomous 
perturbations. The time-varying contributions due to the Moon, the planets and other celestial 
bodies must be included in the list of relevant perturbing effects. On the basis of the present work, 
we are confident that the study of a properly modified version of the equations of motion in the 
form of a forced system of coupled nonlinear oscillators (e.g. that recently presented in ref.|12]) 
will provide additional solutions with a degree of accuracy comparable to that presented here. 
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